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I. INTRODUCTION 



A Hamiltonian system of A degrees of freedom is said to be completely integrable, in the 
Liouville-Arnold sense, if it possesses functionally independent globally defined and single- 
valued A integrals of motion in involution [fj],@. It is called to be superintegrable if it admits 
more than A integrals of motion. Not all the integrals of superintegrable system can be in 
involution, but they must be functionally independent otherwise the extra invariants are 
trivial. In analogy to the classical mechanics, a quantum mechanical system described in 
A-dimensional (AD) Euclidean space by a stationary Hamiltonian operator H is called to 
be completely integrable if there exists a set of A — 1 (together with H, A) algebraically 
independent linear operators Xi,i = 1,2,..., A — 1 commuting with H and among each 
other |5|-|TT|]. If there exist k additional operators Yj,j = 1,2, ...,k where < k < A — 1, 
commuting with H it is said to be superintegrable. The superintegrability is said to be 
minimal if k = 1 and maximal if k — N — 1. 

Classical and quantum mechanical examples of the maximally superintegrable systems 
for any finite A are; the Kepler- Coulomb problem, the harmonic oscillator with rational 
frequency ratio, the Calogero-Moser system in a harmonic well and the Winternitz (or, 
Smorodinsky-Winternitz) system. The first two are known, for A = 3, since the time of 
Laplace and the superintegrability of the last two systems were established for the first time, 
respectively, by Wojciechowski || and by Evans 0. The first systematic search for other 
possible superintegrable systems was begun by Winternitz and coworkers. They firstly found 
four independent 2D potentials that are separable in more than one coordinate system ||, 
and then they extended this to A = 3 0]. This approach is based on two assumptions; 
(1) Hamiltonians are of potential form. (2) Integrals of motion are at most quadratic in 
momenta (or, in derivatives). The Winternitz program has been completed in Ref. |7j where a 
complete list consisting, up to the equivalence of linear transformations, thirteen different 3D 
potentials with four or five independent integrals of motion is given. Winternitz potentials 
have also been considered by different formulations such as path integral formulation M, 
Lagrangian formalism || and evolutionary vector fields formalism [|1C . 

In this paper we report an infinite family of 2D potentials which are not only superinte- 
grable, but at the same time isospectral. We shall give explicit expressions of the potentials, 
their dynamical symmetry generators and the algebra they obey as well as their degenerate 
spectrum and corresponding normalizable states. We achieve this goal by following an alge- 
braic method which is based on and, in fact, is an extension of intertwining operator idea. 
This is closely connected with the supersymmetric (SUSY) methods such as the Darboux 
transformation, and Schrodinger factorization which deal with pairs of Hamiltonians having 
the same energy spectra but different eigenstates [|T2Ln~5[| . It turns out that each member 
of this infinite family is a triplet of potentials one of which is the same for entire family 
and the other two change from member to member. Hence, we have, in fact, two different 
infinite families of superintegrable and isospectral potentials. The fixed potential turns out 
to be one of 2D Winternitz potentials and determines the spectra of both families and the 
other two are intertwined to it by Darboux type transformations. The generators of these 
transformations depend on eigenf unctions of two associated solvable ID problems that result 
from the separation of the Winternitz potential in different coordinates. We should empha- 
size that our approach makes it possible to apply Darboux transformations simultaneously 
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to potential and to its symmetry generators in such a way that superintegrability property 
is preserved. 

Formal aspects of our method together with a brief review of the main points of the 
intertwining operator idea will be given in the next section. Sees. III-V are devoted to 
explicit realization of our method. In Sec. IV we present the most general form of 2D 
integrable and isospectral potentials in the plane polar coordinates. Two subfamilies of 
superintegrable and isospectral potentials and then their general forms are presented in 
Sees. VI and VII. Sec. VIII contains a review of bound states of the associated ID problems 
and the above mentioned Winternitz potential. After investigating the symmetry generators 
and their algebra in Sec. IX, the normalizable states of the generated superintegrable and 
isospectral potentials are given in Sec. X. 



II. MULTIPLE INTERTWINING METHOD 

The object of the intertwining method is to construct a linear differential operator £ 
which intertwines two Hamiltonian operators Hq and Hi such that £H Q = Hi£. Two im- 
portant facts that immediately follow from this relation are; (i) If ip° is an eigenfunction 
of Hq with eigenvalue of E° then ip l = £ip° is an (unnormalized) eigenfunction of Hi with 
the same eigenvalue E° . (ii) When Hq and Hi are self- adjoint (on some common function 
space) £t intertwines in the other direction Hq& = C^Hi and this in turn implies that 
[H ,&£] = = [CC\Hi], where t and [, ] stand for Hermitian conjugation and commuta- 
tor. The first property shows that £ transforms one solvable problem into another, and the 
second one means that two hidden dynamical symmetries of H and Hi are immediately 
constructed in terms of £. These are dimension and form independent general properties of 



this method [|T^,|T^]. In the context of ID systems where £ is taken to be the first order differ- 
ential operator and Hamiltonians are of the potential forms two additional properties arise; 
(i) Every eigenfunction of H (without regard to boundary conditions or normalizability) 
can be used to generate a transformation to a new solvable problem, (ii) A direct connection 



to a SUSY algebra can be established |18| . The first property is a manifestation of the cele- 
brated Darboux transformation and its generalization (Crum transformation). The second 
property enables us to express in a compact algebraic form of the spectral equivalence of 
the intertwined systems. 

Now suppose that there are three self-adjoint Hamiltonian operators H , Hi, H 2 which 
are intertwined as 

£ioH = Hi£ w , £ 2 \Hi = H 2 £ 2 i. (1) 

The subscripts of the intertwining operators are used to distinguish them and to denote the 
intertwined Hamiltonians. Eqs. (1) immediately imply that £ 2 o = £ 2 i£io will intertwine 
H and H 2 as follows, 

£2qHq = H 2 £ 2 q. (2) 
Eqs. (1) and (2) can be unified into the following diagram 

Hq — > Hi 

\ I (3) 
H 2 



3 



which must be understood in the sense described by (1) and (2). 
Adjoint s of (1) and (2) yield 

£21-^2 = HxL 2 xi C\ Q Hi = H C\ , C,\qH 2 = HqCIo- (4) 

That is, the adjoints of the intertwining operators will intertwine in the reverse directions 
and this can be represented by a diagram the same as (3) with reversed directions of arrows. 
Making use of (1-4) it is easy to show that each of H , Hi, H 2 has two dynamical symmetry 
generators respectively given by; 



(5) 



The subscripts of Xj, Yj indicate the Hamiltonians they belong to. Throughout this paper we 
assume that the domains of definition of Hamiltonians and intertwining operators are some 
linear subspaces of a common Hilbert space Ti = L 2 (Q) with the standard sesquilinear inner 
product. L 2 (Q) is the space of all square-integrable functions (and distributions) defined on 
a subspace Q of iVD Euclidean space R N pT9|— pTj| . 

For all N > 2, the diagram (3) implies a triplet of isospectral Hamiltonians such that each 
has two dynamical symmetries. By construction, all the symmetry operators obtained in 
this manner will be factorized, and have even orders depending on the order of intertwining 
operators. They will be of the same order only for Hi. According to the von Neumann 
theorem (see Ref. (T^j, pp.141 and Ref. [p0| , pp.275) £ijC\j ( and C^Cy) are self-adjoint 



X 




Y 


= £20^20 


Xi 


= £io-£io> 


Yi 


= ^21^21 


x 2 


= ^21^21) 


Y 2 


= £20^20 



and nonnegative if Cij are closed with dense domains of definition. Otherwise there may 
exist states in which they have negative expectation values (see Sec. X). If C w and £21 are 
taken to be algebraically independent, the independence of Xi, Yi pairs will be guaranteed 
from the outset. Extensions of these ideas to higher dimensions will be, generically, called 
multiple intertwining method. A simple observation that this work initiated from is that, in 
the particular case of N = 2 the diagram (3) guarantees the superintegrability of the three 
Hamiltonians. In the case of N = 3 such a diagram will imply, provided that symmetry 
generators are commutative, the integrability of the potentials. 

The rest of the paper is devoted to explicit realization of these formal observations for 2D 
systems. Firstly we will determine the most general form of the potentials and the first order 
intertwining operator for two Hamiltonians. We then construct the intertwinings H — > Hi 
and Hi — > H 2 by special forms of the intertwining operator. We end this section by ex- 
plaining our use of the adjective "isospectral". Two Hamiltonians are said to be isospectral 
if they have the same eigenvalue spectrum [|l^,^,^]. In this sense two linearly intertwined 



Hamiltonians are always formally isospectral except the eigenvalues corresponding to the 
kernel of the intertwining operator. Even for these exceptional cases one can construct 
eigenfunctions corresponding to these eigenvalues at least for ID and 2D systems by appeal- 
ing to the Liouville formula and its 2D version [14]]. However, due to physical requirements, 



in the case of the bound states mainly due to normalizability conditions, some eigenvalues 
of one of the partner potentials are to be discarded. For higher dimensional systems also the 
degree of degeneracy of a common eigenvalue may be different (see Sec. X). These will just 
mean that a finite number of eigenvalues are to be disregarded for they are not physically 
admissible. 
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III. INTERTWINING IN TWO DIMENSIONS 



We start by considering a pair of 2D one particle systems characterized by the Hamilto- 
nian operators of potential form, 

H t = -V 2 + Vi, H f = -V 2 + V f , (6) 

where the potentials Vi, Vf (and eigenvalues of Hi, Hf) are expressed in terms of 2m/ 'h 2 and 

v 2 = % + \a r + 1^ 

is the Laplace operator in the plane polar coordinates (r,9). m is the mass of the particle 
and h denotes the Planck constant. Here and hereafter we use the notation d x for partial 
derivative d/dx and the subindexes V, "/" as the shorthands for the "initial" and "final". 
We suppose that the Hamiltonians are intertwined by 

CfiHi = HfCfi (7) 

and propose the ansatz that Cfi is the most general first order linear operator 

Cfi = L + L d = L + Lid r + L 2 d e , (8) 

where Ld = L\d r + L^de will be referred to as the differential part of Cfi. The potentials 
and L , L\, L 2 are some real functions of (r, 6) which are to be determined from consistency 
equations of the intertwining relation (7). 

In view of (6) and (8) the relation (7) explicitly reads as 

[V 2 , L d ] = -[V 2 , L ] + [V, L d ] + PC fi , (9) 

where P — Vf — Vi. The second order derivatives come, together with some first order 
derivatives, only from 

[V 2 , L d \ = (V% + ^L^dr + (V 2 L 2 )d e + 



2{-d e L 1 + d r L 2 )d e d r + 2(<9 r Li)<9 2 + + rd e L 2 )d 2 e , (10) 



and by setting their coefficients to zero we obtain; 

d e L x + r 2 d r L 2 = 0, d r L x = 0, L x + rd e L 2 = 0. 
It is straightforward to show that the general solutions of these equations are 

Li = Asin(# + 0), L 2 = B + -cos(6 + <f)), (11) 

where A, B and <fi are integration constants. Since V 2 Li = —L\/r 2 , and V 2 L 2 = 0, we have 
from (10) [V 2 , L d ] =0. As a result of this the relation (9) simplifies to 

[V 2 , L ] = -LAV - L 2 d e V t + P(L + L d ). (12) 
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By substituting 



[V 2 , L ] = V 2 L + 2(d r L )d r + ^(d e L )d e , 

into (12), and then by equating the coefficients of the first and zeroth powers of derivatives 
we obtain 

2d r L = PL U (13) 
2d 6 L = r 2 PL 2 , (14) 
(-V 2 + P)L = L^M + L 2 d V t . (15) 

These three partial differential equations, the first two of which are linear and the third 
is nonlinear, constitute a reduced form of the consistency conditions for three unknown 
functions L , Vi and Vf. 



IV. GENERAL FORM OF 2D INTEGRABLE ISOSPECTRAL POTENTIALS IN 

POLAR COORDINATES 

Eqs. (11), (13-14) and the compatibility condition d r d e L = ded r L imply that 

2V 2 L = L d P, ZL = 0, ZP = 2BrP, 

where Z = Lxdg — r 2 L 2 d r . From the second and third of these equations (or, from (13) and 
(14)) we have L = f(w), and P = —2A 2 f'{w)/r 2 L\, where / is an arbitrary function of 

A 

w = Bcot(9 + 4>) + 



r sin (6* + 0) 

Prime stands for derivative with respect to the argument and when there is no risk of 
confusion the argument will be suppressed. By combining 2V 2 L = L d P with (15) and 
using the found L and P we obtain an inhomogeneous equation from the general solution 
of which the general form of potentials are found to be 

Vi = h(K) + ^ Vf = h(K) + Y±M. (16) 

K K 

Here h is an arbitrary function of k — [A 2 + B 2 r 2 + 2ABrcos(8 + 0)] 1//2 such that L d h = 
and 

V ± {w) = f 2 {w)T{w 2 + B 2 )f{w). (17) 

Eqs. (16) represent the most general form of 2D integrable and isospectral potentials in 
polar coordinates. 

Let us define the operators 

Ti = cos 9d r — - sin 9d e , T 2 = sin 9d r + - cos 9d e , J = d e , (18) 
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which close in the defining relations of the Euclidean Lie algebra e(2) in two dimensions 

[J,T 1 ]=-T 2 , [J,T 2 ]=T Xl [T u T 3 ] = 0. (19) 

Now Ld can be rewritten as 

L d = A sin 0T\ + A cos 0T 2 + 5 J, (20) 

which shows that the differential part of Cfi is an element of e(2). In terms of the Cartesian 
coordinates x = rcos9,y = rsin^ we have Ti = d x ,T 2 = d y ,J = xd y — yd x and T- = 
—Ti, J' = — J. These relations can also be verified from (18) by noting that (d r Y = — (r _1 + 
d r ), (deY = —dg. Now from (5) and (17) the symmetry generators of Hi and Hf are 

£/iAf« = V- ~ L% £fi£fi = v+ - L% 
where L\ is at most quadratic operator in generators of e(2). 

V. CONSTRUCTION OF THE INTERTWINING OPERATORS 

We shall construct the legs of the diagram (3) by adopting particular forms of (20) as 
the differential parts of £io and £21- In doing that we shall make use of the orbit structure 
of e(2) under the adjoint action of the Euclidean group E{2) in two dimensions 

Under a unitary similarity transformation, generated by 

JJ — e a, J e a 1 T 1 +a 2 T 2 £jf _ jj— 1 _ e -(a 1 T 1 +a 2 T 2 ) e ~a J {21) 

where a^s are real parameters and U^ 1 stands for the inverse of U G E{2), the relation (7) 
transforms into CfiHi = HfCfi, where X = UXW . Since V 2 = T\ + T 2 2 is the Casimir 
invariant of e(2), only Vi, Vf and Cfi will change under this E{2) action. Now suppose that 
Ld is of the form (20). Making use of the well known operator identity 

b 2 

e bK Me -bK = M + b [ K} M ] + [ K} M]] H , 

where b is a constant and K, M are two arbitrary operators, one can easily show that 

L d = BJ + e aoJ [T^A sin (p - a 2 B) + T 2 {Acos(p + a 1 B)] e - a ° J . 

Hence, if B ^ we can take Ld = BJ by choosing a\ = —Acos<p/B,a2 = Asin<j)/B. On 
the other hand, if B = 0, A ^ we get L d = AT\ (or, L d = AT 2 ) for the choice a = 4> (or, 
a = —0). Therefore, under the adjoint action of E{2), e{2) has two orbits represented by 
J and T 2 . Since L d and cL d belong to the same orbit for c ^ 0, we can choose L d = J for 
Cio and L d = T 2 for £ 21 . In such a case the potentials and L Q will be specified up to the 
adjoint action of E{2). 

For the first leg H Q — > of (3) we take A = 0, S = 1 in Eq. (11) and redefine the 
Hamiltonians as Hi = H Q and Hf = Hi. Hence Li = 0, L 2 = 1 and Eqs. (13-14) imply that 
/-o fi(') aiul 
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C 10 = H9) + d e , P = V 1 -V = ^f(9), (22) 

where / is an arbitrary differentiable function of 9. Noting that V 2 L = f"(9)/r 2 we obtain 
from (15) and (22) 

vi = Mr) + m, Vl = h(r) + YM, (23) 

where h is an arbitrary differentiable function of r and 

V ± (9) = f(9)±f'(9). (24) 

As a result the first H Q — > leg of the diagram (3) has been constructed. 

For the second leg we take B — 0, A — 1, fix the form of H\ and denote it as — Hi. 
We then look for Hf = H 2 such that C 2 \H\ = H 2 C 2 \ and £ 2 i = -^0 + sin0Ti + cos0T 2 . In 
that case from Eqs. (13-14) we get L = g(u) and 

£21 = g(u) + sin(# + 0)<9 r + i cos(^ + 0)9 , (25) 
P = V 2 -V 1 = 2g'(u), (26) 

where g is an arbitrary differentiable function of u — rsin(^ + 0). It only remains to solve 
the nonlinear equation (15) which now takes the form 

d u [g\u) - g\u)\ = sin(fl + <f>)h'{r) + ^[cos(^ + ^(9) - 2sin(fl + <j>)V+(9)], (27) 

where we have made use of V 2 Lq = g"(u) and of the second equation of (23). Note that 
we could have chosen = 0, but since it costs almost nothing we keep in our formulae in 
order to see that action of E(2). 

Since it further restricts the three arbitrary functions specifying the potentials, Eq. (27) 
is the main equation which determines the final form of the potentials. As a consistency 
condition the right hand side of Eq. (27) must be only a function of u. Nevertheless this 
requirement provides us with many possibilities for /, g and h, which are investigated in the 
next two sections. Note that for any solutions of Eq. (27) the potentials will be connected 
to each other as follows: 

V = V l -^f'(9), V 2 = V 1 + 2g'(u), V = V 2 - 2[g'(u) + A^]. (28) 



VI. TWO SUBFAMILIES OF POTENTIALS 

We construct the simplest family of potentials by taking, in (24) and (27) h = (Ai/r 2 ) + 
a, V + = —Ax. These lead us to 

f 2 + f = -Ai, g 2 -g' = -\ 2 , (29) 

where a, Ai, A 2 are some arbitrary constants. Then, by Eqs. (23-24) and (26), we obtain 
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V = -Af + AO + a, V 1 = a, V 2 = 2{g 2 + A 2 ) + a. (30) 



r 



The general solution of g 2 — g' = — A 2 is 



2 = < 



A^ 2 tan(A^ 2- u + ai); for A 2 > 0, 



i 



I (-Aa) 1 / 2 tanhK-As) 1 / 2 ^ + /or A 2 < 0, 

where a± is a constant. The solution of f 2 + f' = — Ai can be directly read from the above 
relation after the replacement (g,u, A 2 ) — > (/, —9, Ai). 

An important point is that, by the usual linearization of the Riccati equation, if we 
substitute 

into (29) we arrive at two ID Schrodinger equations 

- ^"(0) = Ai^(0), -*"(«) = A 2 tt(u). (32) 

While the second one can be considered as a free motion, this is not the case for the first 
since < 9 < 2tt. An appealing case is to consider one, or, both of them as infinite 
square-well problem. Normalized eigenfunctions subjected to boundary conditions, say, 
■0(0) = = il)(2ir) and corresponding eigenvalues are 

M0)=^ 1/2 sm(he), Xi, k = ^, k = 1,2, ... (33) 
Hence fk = (k/2) cot(k9/2) and by virtue of Eqs. (22) and (24) we have 

V ^ = ^W + a ' ^ = k r ot(lkO ) + Se . (34, 

To distinguish the resulting potentials, corresponding intertwining operators and the param- 
eter Ai we have labelled them by the quantum number k. The w-problem can be treated in 
a similar way. In any case, the potentials and transformations among them are generated 
by solutions of these two auxiliary ID problems. The existence of V\ = a explicitly shows 
that the member potentials are isospectral to a 2D free motion. As a result we have found 
a five parameter (a, ai, Ai, A 2 , 0) family of 2D potentials that are generated, in a nontrivial 
way, by two ID problems. 

We specify a second subfamily of potentials by taking, in (23-24) and (27) 



h = At + 7r«r 2 + a, (35) 
r z 2 



and V + = — Ai. These lead us to the same equation as in (29) for / and to the Riccati's 
equation 

g 2 -g'- l -au 2 + A 2 = 0, (36) 



for g. By Eqs. (23-24) and (26) the member potentials are found to be 



1 2 , 2(/ 2 + AQ 
V = -ar H 5 h a, 

2 r z 

V 1 = \ar 2 + a, (37) 

\/ 2 = i«r 2 cos 2(^ + 0) +2g 2 (u) + (a + 2A 2 ), 

where g is any solution of (36) and / is any solution of f 2 + f = — X 1 . 
Now the ansatz (31) for g transforms (36) into 

- *"(«) + -au 2 <H(u) = A 2 *(«), (38) 

which is the well known Schrodinger equation for the ID harmonic oscillator. In that case 
the entire family will have 2D isotropic harmonic oscillator spectrum given by the eigenvalues 

e£ ] = hu(e + i), £ = o,i,2,..., (39) 

which are £+1 times degenerate for a given t. For concrete examples we recall the normalized 
eigenfunctions and corresponding eigenvalues of the ID harmonic oscillator: 

* n («) = N n e~ p2u2/2 H n ((3u), E n = ^-\ 2 ^ = %u(n+ 1 -), n = 0,1,2,..., (40) 

where N n is the normalization constant, H n denote the Hermite polynomials and 

( mw 1/2 a 1/4 /3 1/2 

In writing Eqs. (39-41) we have restored 2m/h 2 into our notation in which the dimension of 
f3 is (length)' 1 . Like A 2 , also V 2 ,C 2 \ and the function g must be labelled by the quantum 
number n: 

9n(u) = = /? 2 u - d u HH n {(3u)\, 

V 2 {n) = [3 4 r 2 cos 2(9 + 0) + 2<?» + a + A(3 2 (n + 1), (42) 

4? = 9n(u) + sin(fl + <p)8 r + - cos(# + <f>)d e . 

r 

For the first three Hermite polynomials H (x) = l,Hi(x) = 2x, H 2 (x) = Ax 2 — 2 we have 

g = Pu, 9l =f3u--, g 2 = f3 ^-^—^ u. 
Considering the /-problem as above 4? is given by (34) and Vq is 
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VII. GENERAL FORM OF THE POTENTIALS 



Returning the general discussion of Sec. V, the most general potentials are obtained by 
choosing, in Eq. (27), h as in (35) and by postulating the equation 

cos(0 + 0)Vj(0) - 2 sin(0 + <f>)V + [6) = 2X 1 sin(0 + 0) - * (44) 

for V+. It is not hard to check that (35) and (44) are the most general relations which make 
the right hand side of Eq. (27) only a function of the u variable. The general solution of 
Eq. (44) is 

v ^^ + f^-^4w) + ^iwr Xi ' (45) 

where Xi,b and c are some constants. When (35) and (44) are inserted into (27) we obtain 
a new Riccati's equation for g{u) 

g 2 -g' =\au 2 + ^-X 2 . (46) 
2 u 2 

By virtue of (23), (28), (35), (45) and (46) the corresponding potentials can be written as 

V =- 1 r & | C 1 | 2 (/ 2 + A l) | 

2 aT r 2 W(0 + 0) sin 2 (0 + 0) J + r 2 + °' 

VI = + Xos 2 (0 + 0) + sin 2 (0 + 0) ] + °' (47) 

1 X I) C (X 

V 2 = -ar 2 cos 2(9 + 0) + -[ ^— tt - . 2 * ,J + 2(g 2 + A 2 + ^). 

2 r 2 cos 2 (9 + (p) sin (0 + 0) 2 

Vi is immediately recognized as one of 2D Smorodinsky-Winternitz potentials which 
accepts separation of variables in the Cartesian, polar and elliptic coordinates. Being fixed 
in the whole family it determines the structure of spectrum of all potentials. While Vo is 
separable in the plane polar coordinates V 2 is separable only in the Cartesian coordinates. 
Vq and V 2 represent two families of the superintegrable and isospectral potentials generated 
by the functions / and g which are subjected to Eqs. (45) and (46). The normalized 
eigenfunctions, corresponding eigenvalues and the symmetry generators will be the subject 
of the next three sections. 

Having specified the most general forms of the potentials we now show how to develop 
a hierarchy of the potentials. 

On substituting (31) into (45) and (46) we arrive at the following two ID problems 

HptM&) = KkMQ), Hso^niu) = X 2 , n ^ n (u), (48) 

where k and n are possible quantum numbers and 

HpT = —Ja2 + VpT > VpT = 2(1,^ + • 2777 , 7v ( 49 ) 
d9 2 cos 2 (9 + 4>) sin (0 + 0) 

d 2 1 c 

H SO = + Vso > Vso = 2 au2 + ^ 50 ^ 
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These are the well-known generalized Poschl- Teller (PT) and singular oscillator (hence the 
subscript SO), or the radial oscillator potentials. By virtue of (28) and (31) the potentials 
can be rewritten as 

Vi k) = V x -^d 2 e \n^ k {e), Vt } = V 1 -2%hxHf n (u). (51) 

Here and here after we label the potentials by the quantum numbers of the associated ID 
problems that generate them. (51) explicitly shows that Vq^ and are generated from V\ 
by the Darboux type transformations. The functions that generate these transformations are 
the eigenfunctions of the associated ID problems. This constitutes an extension of Darboux 
transformations for 2D problems. Another point worth emphasizing is that any solution of 
these ID problems can be used in generating the potentials. But, as easily accessible results 
from the literature, only normalizable solutions of these problems will be presented below. 
From now on we take = and in Sees. VIII and X we include 2m/ h 2 into our notation. 



VIII. BOUND STATES OF THE ASSOCIATED PROBLEMS AND Vy 

Provided that c > —1/4, the bound states of Hso belonging to the Hilbert space L 2 (0, oo) 
are given as follows |3|,|26|-p8| 

q?l( u ) = N^°J +£V e-^ u2 / 2 L^((3 2 u 2 ), 

h 2 

K = —\\ n = hu{2n + ev + 1), n = 0, 1, 2, (52) 
N so = [ 1 I 1 / 2 u = -(l + 4c) 1/2 

where is the normalization constant, Lff(z) are the generalized Laguerre polynomials, 
P is defined by Eq. (41), T stands for the Gamma function and e — ±. \E^j(u)'s satisfy the 
orthogonality relation p9fl 

POO 

/ %(u)ty £ n ,(u)du = 5 nn ,, (53) 
Jo 

which is valid for ev > —1. This implies that for cG / = [-1/4,3/4) (that is for — 1/4 < c < 
3/4 ) both values of e — ±, and for c > 3/4 only e = + can be used for each n. Although 
the generated potentials do not depend on the normalization constants of the associated ID 
problems we write them for completeness. 

From the most general point of view and in accordance with the fact that Hso is parity 
invariant, defined parity states of H so belonging to the Hilbert space L 2 {— oo, oo) can be 
given as follows [^ 



f \u\ 1 2+^e-^ u2 l 2 L e ^{f3 2 u 2 )- for u > 0, 
*n(«) = T^N™ (54) 
I -e\u\^ +eu e-^ u2 / 2 L^(p 2 u 2 )- for u < 0. 

These obey the following orthogonality relation 
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/oo 
%(u)%,(u)du = 6 nn ,5 eg , (55) 
-oo 

where e, e may equal ±. For e = e (55) follows from the orthogonality of the generalized 
Laguerre polynomials and for e ^ e from the parity reasons as can be verified directly 
from (54). The corresponding energy eigenvalues are given by (52). For c < —1/4 the 
energy spectrum is not bounded from below which implies "falling of the particle to the 



center" and physical interpretation is lost [27,28]. As c — > 0, v — > 1/2 and ^(u)'s go over 



for e = + to odd parity and for e = — to even parity harmonic oscillator wave functions. 
This follows from the relations between the Hermite and Laguerre polynomials [p7|29 . The 
corresponding limits of the energy eigenvalues are obvious from (52). 

The normalized eigenfunctions and corresponding eigenvalues of V\ can now be written 

as 

Ef = E s ni + E £ n2 = hu(2£ + ev + ev + 1), (56) 
£ = ni + n 2 , £, ni,n 2 = 0, 1,2, ... 

where v = (1 + 46) 1/2 /2. ^(z), \E^ 2 (y) and E E nv E e n2 are given by (52) (or, (54)) with 
suitable replacements of the parameters and variables. It follows that bound states of V\ 
exist for b,c > —1/4. For b,c £ I there are four different states for each value of I. In the 
case of6G/,c>3/4, orcG/, 6>3/4 there are two different states for each value of £, and 
one state in the case of b, c > 3/4. In each case, for a given value of I the state with energy 
Ef 6 is (£ + l)-fold degenerate. We should also note that if we require the wavefunctions 
to be separately continuous at the origin the interval / = [—1/4,3/4) and the conditions 
b,c> 3/4 must be replaced as / = [—1/4, 0] and as b, c > 0. 

The singular oscillator problem has the spectrum generating algebra su(l, 1) = { Jo, J± : 
[Jo, J±] = =t J±, [J+, J-] = — 2J } realized as [p6| , |30|l 



Jo = |f, J ± = -\\nuT^d u f-^l (57) 

with the Casimir invariant C 2 = — J + J_ + Jq — Jo = (4c — 3) / 16. Therefore, as will be shown 
in the next section, the symmetry algebra of the ifi-problem is closely connected with this 
kind two commuting copies of su(l, 1) algebra. 

For later use it will be convenient to consider iJpT-problem in relation with the solution 
of /Jx-problem in the polar coordinates. In this case the eigenvalue equation of Hi separates, 
by taking ^^\r, 9) = Ri Cl (r)ip/ : (9), into the Poschl- Teller problem given by (48) and into 
the radial equation 

[-(^ + ~) + p 2 + ^}RM = XRM, (58) 

where p = (3r and A = E//3 2 . In terms of v — sin 2 9, and ijjk(9) = v^^ +ev \l — v)^^ +eu ^ F{v) } 
the eigenvalue equation of H PT leads us to the hypergeometric equation for F(v): 

d 2 F dF 
< l ~ v ^^¥ + [C - ^(7 + V + 1)1-^7 - IVF = 0. 
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The general solution of this equation is 

F(v) = A 2 F 1 ( r/ , 77; C; v) + Bv'-^F^ - c + 1, v - C + 1; 2 - C; v), 

where A and B are arbitrary constants, 2 F\ denotes the hypergeometric function and 

1 1 

7 = -(1 + ev + ev + y/X 1>k ), V= -(1 + ev + ev- \/Ai,fc), C = 1 + £f- 

For normalizable solutions S must be zero and 7 (or 77) must be a negative integer, say, —k. 
In that case the hypergeometric function goes over to Jacobi polynomials P[f u ' £U \l — 2v) 
and the resulting eigenfunctions and eigenvalues can be written as follows |7|J31 

[9) = N[ T sin^ cos5 +e " p 6Pt' £ ~ P) (cos 29) , 

E k = — Ai,* = ^-(2k + ev + ev + l) 2 , = 0, 1, 2, (59) 

2m 2m 

PT _ , 2(2fc + ev + ev+ l)r(fc + l)r(fc + ei/ + ev + 1) 1/2 
fc _i T(k + ev + l)r(Jfe + ev + 1) J ' 

Substituting Ai^ = (2/c + ez/ + + l) 2 into Eq. (58) and trying the solution R kl (p) = 
p^ e -p / 2 Gk 1 (p) we end up, for p, = y/\i tk = (2k + ev + ev + 1), with the equation 

+ (a* + 1 - z)^I - i[2(Ai + 1) - A]G &1 = 0, (60) 

where z = p 2 . Provided that — [2(/j+ 1) — A]/4 is an integer, say fci = 0, 1, 2, the solutions 
of (60) are the generalized Laguerre polynomials. Hence, the radial solutions are 

R kl {p) = Nk^e-^Ll( P % N kl = [ *l ik * + l\j 1/2 . (61) 

1 [ki + p + 1) 

One can easily verify that ^.(0)'s and R kl (pys obey the following orthogonality relations 

poo /"^/^ 

/ R kl {p)R k[ {p)rdr = 6 klk[ , / ^ k {9)^ w {e)d9 = 6 W . (62) 

As a result the eigenfunctions of Hi can be written in polar coordinates as follows 

^« £ " £ (r, 9) = N kl N^{l3rYe-^' 2 sin^ 9 cos^ +£ " p 0L£ (p 2 r 2 )P^ £P) (cos 20), (63) 

with £ = k\ + ^2; ^, ki, k 2 = 0, 1, 2, and /1 = (2/c 2 + ez/ + £z> + 1). Since ^(0) given by (59) 
will be used in generating Vq potentials, in writing (63) we have changed the quantum 
number k as k 2 . Observe that a similar change [n — > n 2 ) has been made in writing (56). 
Note also that the condition — [2(p + 1) — A]/4 = k\ gives the eigenvalue (56) for the V\- 
problem, with I = k\ + k 2 . We should also note that, as has been done in Eq. (54), the 
solutions (63) may be extended to all xy-plane such that they have definite parity under 2D 
parity transformation: (r, 9) — > (r, 9 + 7r). 

Inserting ip k and \fr n into (51) explicit expressions of the potentials labelled by the quan- 
tum numbers k and n become available. Besides that presented in the Sec. VI several more 
special subfamilies can be identified. In doing that one should take care of the range of the 
parameters and the domain of definition for potentials. The bound states of and V 2 n ^ 
will be taken up in Sec. X after considering the symmetry generators in the next section. 
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IX. SYMMETRY GENERATORS AND THEIR ALGEBRAS 



As is apparent from previous two sections, the intertwining operators, symmetry gener- 
ators, and the Hamiltonians H ,H 2 must be labelled by the quantum numbers (k,n) of the 
associated potentials. In terms of e(2) generators the labelled intertwining operators are 

4? = /*(*) + 4 d$* = f k (o)-j, 

C^=g n (u)+T 2 , 4f = g n {u) - T 2 . (64) 
It is easy to verify that they obey the following commutators 

[c[tc[f] = 2f k (e), [4i ) ,4i )t ] = 2^H, 

[4?, 4?] = K - k,n) + T u [4o f , 4i )f ] = ~K { - ' n) + T u (65) 

[^10 5^21 J — A + — J l> [Mo i L '2l\ — ~- tV + ~ J l> 



where 



xi fe ' n) =rcos%»±^ 
By virtue of (28) we have 

K[ h ' n) = l -rcos6[H { 2 k) - if< n) ], *l fe ' n) = ±r cos 0[if£° + H { 2 n) - 2^]. (66) 
It will be convenient to start with the symmetry generators of H\ 

A l —^10^10 — n PT — M,k, {VI) 
Y l — L 2\ L 2\ — n SO — *2,n, (OoJ 

where we have made use of T 2 g n {u) = g' n (u). H PT and H S o are defined by (49) and (50). 
Second order symmetry generators of and H 2 n ^ can also be written as follows 

X™ = 4o )t 4o ) = - J 2 = SU - Kk, (69) 

= 4" ) 4i )t = gl + g'n- Tl = m - A 2 , n , (70) 

where 

R pt = -§2- + V pt ~ 2^1n^W, (71) 

42 = + V SO - 2d 2 u In (72) 

These are the so called super partners of Hpt and Hso- As a result, the Hamiltonians of 
ID auxiliary problems are, up to some constants, the second order symmetry generators of 
Hi and their super partners are the second order symmetry generators of and H 2 n \ 

The simplest forms of remaining fourth order generators seem to be their factorized forms 
given by (5). Making use of (65) and (67-70) these can be expressed in a variety of ways, 
some of which are as follows; 



1 
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v (k,n) _ r {k)\ v {n) r (k) _ Ak)\ tt Ak) x x (k) 

= Y} n) X ( k) - [(K {k ' n) - 71)4? + 4? f (tf? ,n) + ^i)]4cn (73) 
= V 1 (n) - £f )f [(^' n) - 71)4? + 4? f (^i fc ' n) + 71)], 

v (k,n) _ r (n) v (k) r (n) J f _ r (n) tt r {n)\ x v (n) 

= X {k) X^ - [(K (k ' n) + 71)45" " C { $(Kf n) + 71)]4? f , (74) 

.(*)■ 

'10 



x W x (fc) _ C ^[C { S{K^ n) - 71) - (i^' n) - 71)4 fe)tl 



At this point we have to emphasize the followings. The existence of Ai^ and A2, n as additive 
constants in x[ k \ and X^^X^ seems to be redundant in regard of superintegrability 
of Hi. In particular, our labelling of the fourth order generators with two indices may seem 
as if we have more symmetries than is needed for superintegrability. But an inspection of 
the first lines of Eqs. (73) and (74) immediately shows that, for a given, say, k and all n the 
set {YQ k ' n \ Xq^} spans only a 2D vector space. As fourth order symmetries labelled with 
one index one may take 

V-( fc ) — v( fe ' n ) _l_ \ — AW Tf Ak) 

■Tr{n) _ v (k,n) , v (n) _ An) tt /*(«)t \ 

r 2 = r 2 + A l,fc A 2 — ^21 n PT>~"l\ ~ A l,fc A 2 • 

However, for overall consistency of the hierarchy such as intertwining of H\ with Hq®, 
and, as we will show in the next section, in determining the spectra of and H^ these 
seemingly redundant constants and labels play an essential role. 

One of virtues of our approach is that the commutativity of the symmetry generators 
with the corresponding Hamiltonian is guaranteed by construction from the outset. For 
justification we first note that the relations 

[H u x{ k) ] = 0=[H u Yl% A; = 0,1,..., (75) 

immediately follow from the fact that x[ k ^ and Y^ emerge from the separation of Hi in 
different coordinate systems. Secondly at a glimpse of Eqs. (67-72) we observe that 

x (k) = X (k) _ 2d 2 lnMe)j X (n) = yin) _ 2d 2 Mn{u y 

That is, the second order symmetry generators of H^ and H^ are Darboux transforms 
of symmetry generators of Hi as are, apart from the factor r~ 2 , H^ and H^ Darboux 
transforms of Hi along different legs of diagram (3). In view of this fact the relations 

[tf.tf] = = [H^ n \xt\ n, k = 0, 1, (76) 
follow from, or, in a sense, are Darboux transforms of (75). Only the explicit check of 

[Hg*\ Yf ' n) ] = = [Ht\ Y} k '% n, k = 0, 1, (77) 

takes tediously a lot of time. This shows an advantage of our method compared with 
the conventional approach in which much effort is devoted to verify the commutativity for 
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specified forms of generators. There it is known that for symmetries higher than second 
order, equations resulting from commutativity are almost intractable. 

It is not so hard to check that [X, , Y}°] ^ 0, j = 0, 1, 2, since the highest order derivatives 
with constant coefficients will appear at the right hand side. For example, 

[Xjp, Y (k ' n) ] hot = [J 2 , JT 2 J] = 4TxT 2 J 3 + 2(2T\ 2 - 2T 2 2 + T{T 2 )J 2 - 821 T 2 J, 
[Xi n \ Y 2 {k '% ot = [T|, T 2 J 2 T 2 ] = 2T 2 4 - 6T 2 T 2 2 - AT^J, 

where [XQ k ' n \ Y^'^hot represents only the highest order terms resulting from [XQ k \ Yq*' 7 ^] . 
Therefore, the symmetry generators of each potential do not close in a finite dimensional 
Lie algebra structure. Note that by the Jacobi identity Zj = [Xj,Yj],j = 0,1,2, is also a 
symmetry generator, but it is algebraically dependent to Xj and Yj. 

There is an elegant way of expressing the symmetry algebra of Hi. For this purpose we 
introduce the generators 

X ± = i L(-^±lax 2 + A) > Di = \(l + 2xd x ), (78) 
Y ^^(- d y^y 2 + ^ D 2 = \(l + 2yd y ), (79) 



which obey the Lie algebras 



[X ± ,Di]=X T , [X + ,X_\ = D 1 , (80) 
[Y ± ,D 2 ] = Y T , \Y + ,Y.] = D 2 , (81) 

with the Casimir invariants 

Xl-X 2 +Dl = ^^, Y 2 -Y 2 + D 2 2 = ^^. (82) 
16 ' + " 2 16 K J 

It is straightforward to show that in terms of (78) and (79) we have 

Hi = 4[3 2 (X + + Y + ), 
X[ k) = 8(X + Y + - X_Y- + DiD 2 ) + K, (83) 
Y} n) = 4f3 2 Y + - A 2>B , 

where K = b + c — Ai^ — (1/2). Eqs. (80-82) are defining relations of two commuting copies 
of a su(l, 1) algebra which can be written as a direct sum su(l, 1) © su(l, 1). The basis 
given by (78-79) is connected with that mentioned in Sec. VIII by linear transformations, 
for instance, by comparing (57) and (79) we have Y + = J ,Y_ ± D 2 = J±. Eqs. (83) show 
that the symmetries of Hi are quadratic in the generators of centrally extended (because of 
the constant K) su(l,l)©su(l,l) algebra. 
By defining 

W (k,n) = l[ X { k \ Yi n) ] = 4(3 2 {X_D 2 - Y-Di), (84) 
8 

one can easily show that 
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[ X { k \ W^} = {x[ k \ Y} n) } + x[ k) (2\ 2 , n - H 1 ) + 

(2Y{ n) + 2A 2 , n - H 1 )(Xi ik - 1) + H^b - c), (85) 
[ Y } n \ W^] = (y/ n) + A 2 , n )(i? 1 - - A 2>n ) - 2/3 4 (xf ) - if), (86) 

where {, } represents the anti commutator. These explicitly show that the extended symme- 
try algebra of H\ spanned by {Hi, x[ k \ Yj. , W^ k,n ^}, with the inclusion of \y( k,n \ closes in a 
quadratic associative algebra for all values of k, n. We also observe that this algebra is a cu- 
bic associative algebra in the enveloping algebra of the centrally extended su(l, 1) (Bsu(l, 1). 
Recently such finitely generated associative algebras have attracted a great deal of interest. 
The structure we have obtained coincides, up to some additive constants, with that presented 



in Ref. JJj] for the Winternitz potential V%. In Ref. |30| this structure is constructed as a 
cubic associative algebra in which counterparts of X 1 ,Y 1 are taken to be purely quadratic 
in the generators of su(l, 1) © su(l, 1). We end this section by emphasizing that exploring 
similar algebraic structures for Hq and and connection between them seems to be an 
important problem which deserves to be taken up in another study. 

X. BOUND STATES OF H^ K) , H^ N) AND THEIR DEGENERACIES 

Representing , , and \Er£^ ES given by (52) and (56), in the Dirac notation, respec- 
tively by the kets \n\e >, \n 2 e > and |1; £ee >, we write (56) as follows 

\l;£e£ >— \riie > \n 2 e >, £ = n\ + n 2 . (87) 

In this notation, the corresponding isospectral states of H 2 are 

\2n;£ss>= £ { 2 r i ) \l;£ee> . (88) 

From (53) (or, 55), (56), (68) and (87) one can easily show that 

< 2n; £ee\2n; £ee > =< 1; £ee\Y} n) \l; £ee >, 

= < n 2 e\H so \n 2 e > -A 2 , n , (89) 
= 2Tiuj{n 2 — n), 

where < -|- > represents the usual inner product of 7i — L 2 (R 2 ) and in the third line we 
have included 2m/ h 2 into the notation. Since £ = n\ + n 2 , this implies that as physically 
acceptable states only those with £ > n will survive in the spectrum of H 2 n \ Moreover, the 
degeneracies of the survived states will be shifted to £ — n since the states corresponding to 
n 2 < n can not be normalized. As a result, the normalized states of H 2 1 ^ are as follows 

\2n;£ee >= [2huj{n 2 - n)]~ 1/2 C { 2 f\l; £ee >, (90) 

provided that £ = n\ + n 2 and n 2 > n. 

In a similar way, if we represent ipk 2 an d -Rfci's given by (59) and (61), respectively by 
the kets \k 2 ee > and \kiee >, the states given by (63) can be expressed as 



18 



\l;£ee >=\kiee > \k 2 ee >, £ = k 1 + k 2 . (91) 

In that case the corresponding isospectral states of Hq are \0k;£ee >= C[ k ^\l;£ee > and 
by virtue of (62), (67) and (91) we have 

< OA;; £ee\0k; £ee > =< 1; £ee\x[ k) \l; £ee >, 

= < k 2 ee\H PT \k 2 ee > -Ai, fc , (92) 
2h 2 



(k 2 -k)(k + k 2 + Ev + Ev + 1). 

m 



Hence, the normalized states of are 

|0Jfe; £ee >= [ (k 2 - k)(k + k 2 + eu + ev + l)]~ 1/2 4o V; >, (93) 

provided that £ = k\ + k 2 and k 2 > k. In this case the degeneracy of the state \0k; £ee > is 
£ — k. Explicit functional realizations of the states (91) and (93) can easily be obtained by 
applying £21 , ^io to the wave functions given by (56) and (63). 

XI. CONCLUDING REMARKS 

The method of intertwining is a unified approach widely used in various fields of physics 
and mathematics such as in investigating particle propagation on a curved space |l6l|2l||33"ll , 
in constructing matrix-Hamiltonian to realize higher dimensional superalgebras [O-pJ 
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solving both ordinary and partial differential equations [16|, in generating exact solutions 



of non-stationary Schrodinger equation [0,RH], and in constructing isospectral potentials 



in an arbitrary space dimension p5| . The method we have introduced increases the power 
and enlarges the range of applicability of the intertwining operator idea. It allows us to 
perform Darboux transformations in higher dimensions in such a manner that, in addition 
to their isospectral deformation property they acquire integrability and superintegrability 
preserving property. In particular, as we have shown the realization of this method for 2D 
systems generates two infinite families of isospectral and superintegrable quantum systems 
intertwined to a 2D Winternitz system. Work on 3D realization of the method is in progress. 

The space of purely second order operators quadratic in the generators of e(2) has, 
under the adjoint action of £"(2), only four orbits whose representatives can be taken to 
be; , J 2 , J 2 + a T 2 and T\ J + JT 1; where a is a constant. Existence of only four types 
Winternitz potentials is closely connected with this orbit structure since each corresponds 
to a different 2D orthogonal coordinate system |3,23]. T 2 , J 2 constitute the differential parts 
of the symmetry generators of V\ and account for its separation in the Cartesian and polar 
(hence in elliptic) coordinates. Therefore, the appearance of the Winternitz potential Vi 
as the common member of two families is of no surprise; it is a direct result of our orbit 
prescription in constructing the intertwinings in Sec. V. We also observe that since only T 2 
and J 2 can be factorized as LC) (or, as C) C) the other three Winternitz potentials can not 
be utilized as V± in the context of this paper. In this regard, a combination of our method 
and the conventional approach may be used for similar purposes. Finally we point out that 
what made it possible to implement Darboux transformations in our approach is that when 
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the eigenvalue equation of V\ is separated in the Cartesian and polar coordinates, at least 
one of the separated equation is of the Schrodinger type. 
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